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EXPERIMENTAL MEASUREMENT AND THEORETICAL MODELING OF 
MICROWAVE SCATTERING AND THE STRUCTURE OF THE SEA SURFACE 
INFLUENCING RADAR OBSERVATIONS FROM SPACE 


INTRODUCTION 


The electromagnetic bias is an error present in radar altimetry of the ocean due to 
the non-uniform reflection from wave troughs and crests. The electromagnetic bias was 
first measured by Yaplee et al. [1970] from an ocean platform using a one nanosecond pulse 
X-band radar. Examples of the normalized radar cross section versus wave displacement 
for a calm and wind driven sea were reported. It was demonstrated that the reflectivity 
in these two cases was not uniform but increased toward the trough. The mean reflecting 
surface was 5% of the significant wave height (SWH) lower than the mean sea surface 
for both examples. The difference between the mean reflecting surface and the mean sea 
surface has been named the electromagnetic bias. 

A study of the electromagnetic bias became necessary to permit error reduction 
in mean sea level measurements of satellite radar altimeters. If not corrected the elec- 
tromagnetic bias could introduce errors in mean sea level measurements as large as 50 
cm [Slinn, 1990]. During 1980 three airborne electromagnetic bias experiments were per- 
formed. Walsh et al. [1984] measured the electromagnetic bias at 36 GHz as 1.1% of the 
gWH. Choy et al. [1984] measured the electromagnetic bias at 10 GHz as 3-5% of the 
SWH. During 1988 Melville et al. [1991] measured the electromagnetic bias at 14 GHz 
as 3.3% of the SWH from an ocean platform. At optical frequencies the electromagnetic 
bias was measured by Hoge et al. [19S4] as biased toward the crests by 2% of the SWH 
for a low wind speed case and biased toward the trough by 0.75% of the SWH for a high 
wind case. Walsh et al. [1989] report additional measurements of the electromagnetic bias 
at optical frequencies for high wind conditions. They found the electromagnetic bias at 
optical frequencies to be unbiased or biased toward the crests as much as 0.5% of the SWH. 
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Satellite radar altimeter measurements have been used to find upper and lower 
bounds for the electromagnetic bias. Studies of the GEOS-3, SEASAT and GEOSAT 
altimeter data [Lipa and Barrick, 1981; Born et al M 1982; Hayne and Hancock, 1982; 
Douglas and Agreen, 1983; Nerem and Shum, 1990; Fu and Glazman, 1991] lead to an 
electromagnetic bias that is about 2-4% of the SWH. 

Jackson [1979] was the first to study the electromagnetic bias theoretically. He based 
his work upon specular point theory [Barrick, 1968] and a model for the joint height-slope 
probability density function based on the work of Longuet-Higgins [1963]. The later work 
of Huang [1984], Barrick and Lipa [1985] and Srokosz [1986] continued the use of specular 
point theory while improving the model used for the joint height-slope probability density 

function. 

During the EM bias experiments by Melville et al. [1990,1991] a wire wave gauge 
was used to obtain the modulation of the high frequency waves by the low frequency 
waves. It became apparent that the EM bias was primarily caused by the modulation 
of the short waves [Arnold et al., 1990]. This report will present a theory using physical 
optics scattering and an empirical model of the short wave modulation to predict the EM 
bias. The predicted EM bias will be compared to measurements at C and Ku bands. 

PREDICTION OF EM BIAS 

The back scattered power from a small patch on the ocean surface depends on the 
displacement of the patch under observation from mean sea level. It has been observed 
that more power is reflected from the troughs of waves than from their crests. A typical 
measurement of the relative back scatter coefficient as a function of displacement is shown 

in figure 1. 

The EM bias can be defined mathemetically as the ratio of the first two moments 

of the back scatter coefficient profile given by 

_ E[r?g 0 (rj)] (!) 
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The task then is to develope a theory to predict the back scatter coefficient profile from 
which the EM bias can be calculated. 


Physical optics scattering will be used to predict the back scatter coefficient. The 
physical optics integral will be evaluated at the frequency of interest rather than the high 
frequency limit. This will have the potential of predicting the electromagnetic frequency 
dependence of the EM bias. 


This theory will assume the primary cause of the EM bias at C and Ku bands to 
be the modulation of the short wave amplitudes by the long waves. This effect will be 
described by a semi-empirical model. The wave number modulation of the short waves 
will be neglected allowing the short wave spectrum to be described by a constant spectral 
shape. The short wave spectrum will be assummed to be unidirectional with a k shape. 
The short wave amplitude modulation will be found empirically by measuring the energy 
in the short waves as a function of long wave displacement by using a wire wave gauge. 


PHYSICAL OPTICS SCATTERING THEORY 


Physical optics or the Kirchhoff approximation is a well known scattering theory 
having been used for rough surface scattering by Beckman and Spizzichino [1963], Hagfors 
[1966], Fung and Moore [1966] and Holliday et al. [1986]. Barrick [1970] argued that the 
only valid use of physical optics was in the high frequency limit resulting in geometric 
optics or specular point theory. Specular point theory was discussed by Kodis [1966], 
Barrick [1968], Tyler [1976] and Barrick and Bahar [1981]. Fung and Chan [1971], Fung 
and Eom [1981] and Chen and Fung [1988] showed using method of moment calculations 
that the high frequency restriction of Barrick [1970] is too restrictive. The validity of the 
Kirchhoff approximation requires only the radius of curvature of the surface to be large 
compared to the electromagnetic wa\elength. 
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The high frequency portion of the ocean wave spectrum causes the average radius 
of curvature of the ocean surface to be small. This apparently renders invalid the use of 
either specular point or physical optics scattering theory. However Tyler [1976] showed 
that the high frequency features of a surface should be smoothed prior to application of 
specular point theory. The common practice is to include only the portion of the ocean 
surface with wavelengths longer than the electromagnetic wavelength [Valenzuela, 1978, 
Barrick and Lipa, 1985]. By numerically evaluating the physical optics integral a filter 
function as proposed by Tyler [1976] does not need to be applied for the employed power 
law spectrum. Physical optics scattering estimates will be compared to method of moment 
calculations for typical conditions of the employed power law spectrum. Physical optics 
scattering estimates obtained without the use of a filter function will be shown to be 
acceptable thus establishing the validity of the scattering theory. 


The physical optics integral for a unidirectional surface is given by (for developement 
see Beckman and Spizzichinio [1963] or Tsang et al. [1985]) 

= (ilL.'j f 1 du( 1 - M) '-^ll-Cp^vLu)) (2) 


where cr° is the back scatter coefficient, k is the electromagnetic wavenumber, L is the 
illumination length, a is the rms wave height, fc( ) is the dominant wave length of the 
surface spectrum and C p is the surface correlation coefficient given by 


fOO 1 

C p (k 0 Lu) = / dk(p-l) A-g L k Pcos(kLu) 


This transform can be performed giving a series solution as 
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The physical optics integral can be evaluated numerically. Figure 2 shows the back 
scatter coefficient as a function of ok for k 0 L = 2tt and p=3. Figure 2 also compares the 
physical optics results with exact method of moment results. Good agreement is achieved 
over the region tested demonstrating the validity of the scattering theory at C band. 
Method of moment code will be improved in the future to allow surfaces with larger ok to 
be tested. 


GULF OF MEXICO EXPERIMENT 

An experiment to measure the EM bias at C and A u bands (the frequencies of 
the Topex/Poseidon altimeters) was conducted from December 1989 through May 1990 
[Melville et al, 1990] from a Shell Offshore production complex (Brazos-19) in 40 meters of 
water off the coast of Texas in the Gulf of Mexico. Nadir looking coherent scatterometers 
at 5 and 14 GHz and a Thorn/EMI IR wave gauge were mounted 18 meters above sea 
level in the middle of a 60 meter bridge joining two platforms. For short periods of the 
experiment, a capacitance wire wave gauge was mounted adjacent to the footprints of the 
scatterometers. The wind speed and direction, air and sea temperature, humidity and rain 
fall were measured by an R. M. Young instrument package. The data contained in this 
report comes from a week of data taken from May 10 through May 17, 1990. 

The EM bias was measured using the back scatter and doppler of the C and Ku band 
scatterometers. The wave displacement was obtained by integrating the doppler centroid, 
which is proportional to the vertical wave velocity, over time to give the displacement. 
This method of obtaining the wave displacement was compared with the measurements 
of the Thom /EMI IR wave gauge and the wire wave gauge yielding close agreement. 
The simultaneous measurements of back scatter and wave displacement were then used to 
calculate the EM bias. 
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The capacitance wire wave gauge was used to measure the short wave modulation. 
The short wave RMS height was measured by calculating the energy in the high pass filtered 
wave gauge output. The wave gauge output was high pass filtered at 1.2 Hz corresponding 
to a one meter separation wave length (assuming a dispersion relation of w 2 = gk.) The 
short wave modulation can then be described by the short wave RMS height as a function 
of displacement from mean sea level. Physical optics scattering and the measured short 
wave RMS height profile can then be used to predict the back scatter coefficient profile 
from which the EM bias can be calculated. 


EM BIAS AND POWER PREDICTION RESULTS 


The wave displacement as measured by the capacitance wire wave gauge is shown 
at the top of Figure 3 for a typical record. This record was recorded on day 6 or May 16, 
1990 at 09:58. The envelope of the high passed filtered wave displacement (short waves) 
is shown at the bottom of Figure 3. It is easily observed that the short waves are being 
modulated by the long waves with the peak modulation occuring near the crests of the 

long waves. 

The short wave modulation can be described by the RMS short wave height as a 
function of wave displacement. This was calculated for a ten minute record, the wave 
displacement shown in Figure 3 was part of this record, and is shown in Figure 4. This 
clearly shows that the short waves have larger amplitudes on the crests of the waves than 
in their troughs. 

Physical optics scattering and the measured short wave RMS height profile can be 
used to predict the back scatter coefficient profile at C and Au bands. Figure 5 shows the 
predicted and measured back scatter coefficient profiles at C and A u bands. The predicted 
back scatter coefficient profiles can be seen to be in close agreement with the measured 
profiles. 
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Hourly averages of significant wave height, wind speed, C and Ku band EM bias 
are shown in Figure 6 for an interval of one week. The predicted bias was calculated 
from the predicted back scatter coefficient profiles. There is close agreement between the 
measured and predicted EM biases. A comparison between the measured and predicted 
biases is shown in Figure 7. The prediction was accurate to within i 2cm over the range 
of observed bias. 

A comparison between the C and Ku band biases is shown in Figure 8. The C 
and Ku band biases Eire the same at small biases but the C band bias becomes larger 
than the Ku band bias as the bias increases. The predicted bias accurately predicts the 
electromagnetic frequency dependence of the measured EM bias at C and A u bands. 

The measured back scatter coefficients are compared to the predicted coefficients. 
Figure 9 shows measured back scatter coefficients at C and Ku bands plotted versus 
the measured rms short wave height. The measurements are compared to the predicted 
relationship given by the solid line. 

The assumption of a unidirectional surface causes the predicted coefficients to be 
higher than the measured coefficients. A constant gain was applied to the measurements 
to account for this difference. A constant gain difference in the back scatter coefficient has 
no effect on the bias. The bias only depends on the relative relationship between the back 
scatter coefficient and the rms short wave height. 

As seen in figure 9 the measured relationship between the back scatter coefficients 
and the rms short wave heights is well predicted by physical optics scattering. This provides 
strong experimental evidence to the validity of the employed physical optics scattering 
theory. 


EM BIAS MODEL 

A model relating the EM bias to the important parameters describing the ocean 
surface will be developed. It will be shown that the bias can be described by a simple 
relationship between wave height, a short wave modulation strength parameter and a term 
depending on the electromagnetic frequency. 
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The short wave modulation profile is approximated linearly by 

^ 2 (»?) = (^n) 2 (l + m_ fc) 

where <Tm is the rms short wave height and m is a measure of the strength of the modu- 
lation. The radax cross section profile can be computed using the physical optics integral 


giving 


= r 1 du(1 _ |u|) e -4{c m kf ll+mn/yf? 

7T 7-1 


Approximating the res profile by the linear relation 




the bias can be determine from 

c = \f? da°(ri) / 10 j 

Jji *°0?) d n *f=° 

The preceding expression can be computed giving the bias in terms of a frequency 
dependent term a , a short wave modulation strength parameter m and the wave height as 

e = -am^ (H) 

where 

_ fn Ml - n)4 (%1) 2 [1 - Cp (kpLu)] e - 4 ( g m^) 2 [ l-C^fepIu)] 

Jq 1 du( 1 - u) e — 4 ( cr m^) 2 [ 1 — Cp(*oi«)] 

The frequency dependent term a is shown in figure 10 for k$L = 2n and p=3.0. 
Data from figure x give ok between 1 and 2 for C band and between 3 and 6 for Ku band. 
Thus the frequency dependent term a is higher at C band than at Ku band causing the 
C band bias to be larger than the Ku band bias. 
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The relationships given by the analytical result of equation (11) are very important. 
This result presents for the first time the specific dependence of the EM bias on ocean 
surface parameters and the electromagnetic frequency. The linear dependence of bias 
on wave height as predicted by this relationship has been well known from experimental 
observations for some time [Walsh et al., 1989]. A correspondence between the short wave 
modulation profile and EM bias was established by Arnold et al. [1990] and reviewed 
in the last section. Equation (11) specifies the correspondence by showing the bias to 
be proportioned to the short wave modulation strength with a proportionality constant 
depending on electromagnetic frequency. 

The data from the one week experiment as described earlier can be used to test 
the theory. Figure 11 shows the a parameter, modulation strength m, and wave height 
normalized C and 7v it band EM biases. 

The a paramter for C band exhibits large fluctuations when the short wave height is 
small. This is due to the ommission of the long waves in the model spectrum. By including 
the long waves in the model spectrum the transition region between small and large short 
wave height as shown in figure 2 becomes smoother causing the parameter to also become 
smoother. The effect of the long waves can be incorporated in the model and this will be 

done in future work. 

The short wave modulation strength parameter m should exhibit the same trends 
as the normalized EM biases. An inspection of the time series in figure 11 shows this to be 
the case. The modulation has more variablity than the biases but it clearly has the same 

trends. 

The wave height, wind speed, measured biases, and biases predicted using equation 
(11) are shown in figure 12. The results are similar to those of figure 6 of the last section, 
but with more variablity. The increase of variablity is due to the linear approximations 
made for the short wave modulation profile and the relative RCS profile. Figure 13 show 
a comparison of the measured to predicted bias. 
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SUMMARY 


Direct measurements of the ocean surface were used in conjunction with physical 
optics scattering to predict the measured EM bias. The predicted bias was accurate to 
within ± 2cm over the range of observed bias. This establishes the short wave modulation 
as the predominant cause of the EM bias at C and A u bands. 

The EM bias has experimentally been found to depend on wave height, wind speed 
and electromagnetic frequency [Melville et ah, 1990, 1991]. The developed EM bias model 
correctly predicted the dependence on wave height, wind speed via the short wave modu- 
lation strength parameter, and the electromagnetic frequency via the a parameter. 

REFERENCES 

Arnold, D.V., W.K. Melville, and J.A. Kong, Theoretical prediction of EM bias, 
Conference Proc. Oceans 90, 253-256, Washington, D.C., 1990. 

Banner, M.L., I.S.F. Jones, and J.C. Trinder, Wavenumber spectra of short gravity 
waves, J. Fluid Mech., 19S, 321-344,1989. 

Barrick, D.E., Rough surface scattering based on the specular point theory, IEEE 
Trans. Ant. Propag., 16, 449-454, 1968. 

Barrick, D.E., Unacceptable height correlation coefficients and the quasi-specular 
component in rough surface scattering, Radio Science, 5, 647-654, 1970. 

Barrick, D.E. and E. Bahar, Rough surface scattering using specular point theory, 
IEEE Trans. Ant. Propag., 29, 798-800, 1981. 

Barrick, D.E. and B.J. Lipa, Analysis and interpretation of altimeter sea echo, Ad- 
vances in Geophysics, 27, 61-100, 1985. 

Beckmann, P. and A. Spizzichino, The scattering of electromagnetic waves from rough 
surfaces, Macmillan Co., New York, N.Y., 503 pp., 1963. 


11 



Born, G.H., M.A. Richards and G.W. Rosborough, An empirical determination of the 
effects of sea state bias on SEASAT altimetry, J. Geophys. Res., 87, 3221-3226, 1982. 

Chen, M.F. and A.K. Fang, A numerical study of the regions of validity of the Kirch- 
hoff and the small-perturbation rough surface scattering models, Radio Science, 23, 
163-170, 1988. 

Choy, L.W., D.L. Hammond and E.A. Uliana, Electromagnetic bias of 10 GHz radar 
altimeter measurements of MSL, Mar. Geod., 8(1-4), 297-312, 1984. 

Desanto, J.A. and G.S. Brown, Analytical techniques for multiple scattering from 
rough surfaces, Progress in Optics, 23, 1-62, 1986. 

Donelan, M.A., J. Hamilton and W.H. Hui, Directional spectra of wind-generated 
waves, Phil. Trans. R. Soc. Lond. A, 315, 509-562, 1985. 

Donelan, M.A. and W.J. Pierson, Jr., Radar scattering and equilibrium ranges in 
wind- generated waves with application to scatterometry, J. Geophys. Res., 92, 49 <1- 
5029, 1987. 

Douglas, B.C. and R.W. Agreen, The sea state correction for GEOS 3 and SEASAT 
satellite altimeter data, J. Geophys. Res., 88, 1655-1661, 1983. 

Evans, D.D. and O.H. Shemdin, An investigation of the modulation of capillary and 
short gravity waves in the open ocean, J. Geophys. Res., 85, 5019-5024, 1980. 

Fu, L. and R. Glazman, The effect of the degree of wave development on the sea state 
bias in radar altimetry measurement, J. Geophys. Res., 96, 829-834, 1991. 

Fung, A.K. and R.K. Moore, The correlation function in Kirchhoff’s method of so- 
lution of scattering of waves from statistically rough surfaces, J. Geophys. Res., 71, 
2939- 2943, 1966. 

Fung, A.K. and H.L. Chan, On the integral for backscattering from a randomly rough 
surface, Proc. IEEE, 59, 1280-1281, 1971. 

Fung. A.K. and H.J. Eom, Note on the Kirchhoff rough surface solution in backscat- 
tering, Radio Science, 16, 299-302, 1981. 


12 



Hagfors, T., Relationship of geometric optics and autocorrelation approaches to the 
analysis of lunar and planetary radar, J. Geophys. Res., 71, 379-382, 1966. 

Hayne, G.S. and D.W. Hancock, III, Sea-state-related altitude errors in the SEASAT 
radar altimeter, J. Geophys. Res., 87, 3227-3231, 1982. 

Hoge, F.E., W.B. Krabill and R.N. Swift, The reflection of airborne UV laser pulses 
from the ocean, Mar. Geod., 8(1-4), 313-344, 1984. 

Holliday, D., G. St-eyr and N.E. Woods, A radar ocean imaging model for small to 
moderate incidence angles, Int. J. Remote Sensing, 7, 1809-1834, 1986. 

Huang, N.E., S.R. Long, L.F. Bliven and C. Tung, The non- gaussian joint probability 
density function of slope and elevation for a nonlinear gravity wave field, J. Geophys. 
Res., 89, 19G1-1972, 1984. 

Jackson, F.C., The reflection of impulses from a nonlinear random sea, J. Geophys. 
Res., 84, 4939-4943, 1979. 

Jahne. B. and K.S. Riemer, Two-dimensional wave number spectra of small-scale 
water surface waves, J. Geophys. Res., 95, 11531-11546, 1990. 

Kodis, R.D., A note on the theory of scattering from an irregular surface, IEEE Trans. 
Ant. Propag., 14, 77-82, 1966. 

Lipa, B.J. and D.E. Barrick, Ocean surface height-slope probability density function 
from SEASAT altimeter echo, J. Geophys. Res., 86, 10921-10930, 1981. 

Longuet-Higgins, M.S., The effect of non-linearities on statistical distributions in the 
theory of sea waves, J. Fluid Mech., 17, 459-480, 1963. 

Melville, W.K., D.V. Arnold, R.H. Stewart, W.C. Keller, J.A. Kong, A.T. Jessup and 
E. Lamarre, Measurements of EM bias at Ku and C bands, Conf. Proc. Oceans 90, 
1S1-186, Washington, D.C., 1990. 

Melville, W.K., R.H. Stewart, W.C. Keller, J.A. Kong, D.V. Arnold, A.T. Jessup, 
M.R. Loewen and A.M. Slinn, Measurements of electromagnetic bias in radar altime- 
try, J. Geophys. Res., 96, 4915-4924, 1991. 


13 



Nerem, R.S., B.D. Tapley and C.K. Shum, Determination of the ocean circulation 
using GEOSAT altimetry, J. Geophys. Res., 95, 3163-3179, 1990. 

Reece, A.M., Jr., Modulation of short waves by long waves, Boundary Layer Meteorol., 
13, 203-214, 1978. 

Rodriguez, E., Y. Kim and J.M. Martin, The effect of small- wave modulation on the 
electromagnetic bias, J. Geophys. Res., 97, 2379-2389, 1992. 

Shemdin, O.H., H.M. Tran and S.C. Wu, Directional measurement of short ocean 
waves with stereophotography, J. Geophys. Res., 93, 13891-13901, 1988. 

Slinn, A., M.S. Thesis, 1990. 

Srokosz, M.A., On the joint distribution of surface elevation and slopes for a nonlinear 
random sea with application to radar altimetry, J. Geophys. Res., 91, 995-1006, 1986. 

Tyler, G.L., Wavelength dependence in radio- wave scattering and specular-point the- 
ory, Radio Science, 11, 83-91, 1976. 

Valenzuela, G.R., Theories for the interaction of electromagnetic and ocean waves - 
a review. Boundary Layer Meteorol., 13, 61-85, 1978. 

Walsh, E.J., D.W. Hancock, D.E. Hines and J.E. Kenney, Electromagnetic bias of 36 
GHz radar altimeter measurements of MSL, Mar. Geod., 8, 265-296, 1984. 

Walsh, E.J., F.C. Jackson, E.A. Uliana and R.N. Swift, Observations on electro- 
magnetic bias in radar altimeter sea surface measurements, J. Geophys. Res., 94, 
14575-14584, 1989. 

Walsh, E.J., F.C. Jackson, D.E. Hines, C. Piazza, L.G. Hevizi, D.J. McLaughlin, 
R E. McIntosh, R.N. Swift, J.F. Scott, J.K. Yungel and E.B. Frederick, Frequency 
dependence of electromagnetic bias in radar altimeter sea surface range measurements, 
J. Geophys. Res., 96, 20571-20583, 1991. 

Yaplee, B.S., A. Shapiro, D.L. Hammond, B.D. Au and E.A. Uliana, Nanosecond 
radar observations of the ocean surface from a stable platform, IEEE Trans. Geosci. 
Electron., 9, 170-174, 1971. 


14 



FIGURE CAPTIONS FIGURE CAPTIONS 


1. Typical measured Ku band relative back scatter coefficient profile. 

2. Physical optics and method of moments back scatter coefficient for k$L = 27r and p 
= 3 versus ok. 

3. Modulation of short waves. 

4. Short wave RMS height profile. 

5. Measured and predicted back scatter coefficient profiles at C and A u bands. 

6. Hourly averages of significant wave height, wind speed, measured and predicted C 
and Ku band EM bias. 

7. Comparison of measured and predicted EM bias. 

8. Comparison of C and A u band EM bias. 

9. Measured and predicted back scatter coefficients versus ok. 

10. Frequency dependent term a for k$L = 2ir and p = 3. 

11. Hourly averages of a parameter, modulation strength m, and normalized C and Ku 
band biases. 

12. Hourly averages of significant wave height, wind speed, measured biases and model 
predictions. 

13. Comparison of measured bias and model predictions. 
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Ku Band Relative Back Scatter Coefficent 



Figure 1. Typical measured Ku band relative back scatter coefficient profile. 
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Figure 2. Physical optics and method of moments back scatter coefficient for IcqL 
2tt and p = 3 versus ak. 
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Figure 3. Modulation of short waves. 
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Figure 4. Short wave RMS height profile. 
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C Band Relative Bock Scatter Coefficient K u Bond Relative Back Scatter Coefficent 


Figure 5. 




Measured and predicted back scatter coefficient profiles at C and K u bands. 
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Figure 7. Comparison of measured and predicted EM bias 
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Figure 10. Frequency dependent term a for IcqL — 2tt and p — 3. 
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Figure 11. Hourly averages of a parameter, modulation strength m, and normalized 
C and Ku band biases. 
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Figure 12. Hourly averages of significant wave height, wind speed, measured biases 
and model predictions. 
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Figure 13. Comparison of measured bias and model predictions. 
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